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METHODS AND APPARATUS FOR MEASURING 
NAVIGATIONAL PARAMETERS OF A 
LOCOMOTIVE 

CROSS REFERENCE TO RELATED APPLICATIONS 

This application claims the benefit of U.S. Provisional Application No. 
60/153,304, filed September 10, 1999, which is hereby incorporated by reference in 
its entirety. 

BACKGROUND OF THE INVENTION 

This invention relates generally to methods and apparatus for 
measuring navigation parameters of a locomotive, and more specifically to methods 
of determining movement and direction of a track-bound transportation apparatus 
using GPS satellites. 

Modem algorithms for train navigation incorporate very sensitive turn 
rate gyroscopes (gyros) that permit detection of low speed .through high number 
turnouts such as a #20 equilateral turnout. These gyroscopes are quite expensive, and 
suffer from coning problems due to superelevation. In addition, measurement drift 
accumulates over time. This drift must be corrected periodically by reference to data 
in a track data base. It would therefore be desirable to replace gyroscopes (gyros) in 
such applications altogether or to provide an alternate system for determining heading 
accurately for low speed and high speed turnouts. Where trains move over a curving 
track, the rails wear down due to friction between the wheels of the train and the rails. 
As used herein, a curve is represented by any change in path direction. Therefore, 
curving means that a path direction is changing. One solution to this problem is to 
provide lubricant to the rails or to the wheels when a train is moving. However, this 
solution is wasteful of lubricant. It would therefore be desirable to provide a method 
and apparatus to more efficiently reduce rail wear on curving tracks by determining 
when a train is on a curve and by applying lubricant. 
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BRIEF SUMMARY OF THE INVENTION 

There is therefore provided, in one embodiment, a method for 
determining motion and/or location parameters of a locomotive. This method 
embodiment includes steps of phase locking a pair of receivers having antennas at 
different locations on a locomotive to reference signals received from a set of at least 
three satellites, determining a set of phase differences between the reference signals 
received at the antennas, and determining at least one of an accurate heading, and/or 
heading rate of the locomotive using the set of phase differences between the 
reference signals. 

The above embodiment can be implemented on a train using 
inexpensive satellite receivers and an inexpensive processor, and can be used in 
conjunction with a gyroscope having reduced specifications to provide accurate 
heading information in both low and high speed turnouts. In addition, the heading 
information provided by this embodiment can advantageously be used to reduce usage 
of lubricants applied to reduce track wear. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 is a drawing of a locomotive showing a relative location of 
receiving antennas of one embodiment of the present invention installed thereon; 

Figure 2 is a diagram of an embodiment of a heading and attitude 
system of the present invention; 

Figure 3 is a top view of a locomotive on a curving track; 

Figure 4 is a drawing of a locomotive engine showing an embodiment 
of a lubricant dispenser applying a lubricant to track rails; 

Figure 5 is a block diagram of an embodiment of a rail lubricant 
control system of the present invention employing a gyroscope and a tachometer; 



WO 01/18558 



PCT/US00/24807 



Figure 6 is a block diagram of an embodiment of a rail lubricant 
control system of the present invention employing two GPS satellite receivers and a 
tachometer; and 

Figure 7 is a block diagram of an embodiment of a rail lubricant 
control system of the present invention employing a lateral accelerometer and a 
tachometer. 

DETAILED DESCRIPTION OF THE INVENTION 

In one embodiment and referring to Figure 1, inexpensive carrier phase 
satellite, for example Global Positioning Satellite (GPS), receivers are used to 
determine train heading very accurately for low and high speed turnouts, which may 
be used in locomotive nav-systems. 

Two antennas 10, 12 are mounted along a longitudinal or lateral axis x 
of a locomotive 14 or train car at a known distance d apart. Antenna 10 and antenna 
12 feed separate carrier phase GPS receiver 16 and GPS receiver 18, respectively. 
Receivers 16 and 18 lock on to reference signals from GPS satellites (e.g., 20, 22), 
and then each receiver 16, 18 performs a carrier phase measurement on a satellite 
reference signal transmission that it is tracking. This measurement consists of a 
fractional phase part <p and an integer count N of phase cycles from an initial epoch t 0 
to epoch t. 

In the following paragraph, a single difference technique is described, 
although the embodiments described herein can also be applied to a double difference 
technique. Referring to Figure 1, antenna 10 of receiver 16 is a distance d away from 
antenna 12 of receiver 18. Both receivers 16 and 18 receive signals from a satellite, 
for example 20, at an elevation 0, which is essentially identical for both receivers 16 
and 18 because of the distance of satellite 20. However, the distance from each 
receiver 16 and 18 from satellite 20 is different by a measurable amount. Letting $ 
denote a phase measurement made by receiver 16 and $ denote a phase 
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measurement made by receiver 18 for a qth satellite (e.g., q=l for satellite 20, q=2 for 
satellite 22), a scaled difference Atf is given by: 

A f = (ft -<t>l) = MNl - M) + c( At! - A/*) + d» h 
where A^ 9 is the single difference for the 4th satellite; 

h is a unit vector in the direction of satellite q, 

d is a vector distance between the two antennas 10 and 12, 



10 



rf> is an integer phase count for the pth receiver (e.g., p=\ or 2 for 
receiver 16 or 18, respectively) and for the <yth satellite; 

Af*, is an error associated with the pth receiver for the <yth satellite; 
c is the speed of light; and 

X is a wavelength of the radio signal received from the qth satellite. 



To solve for d , it is assumed that the clock error and the integer 
ambiguity between both receivers to the satellites are known. The measurements are 
made to each of n satellites in the GPS system, and integer phase counts and clock 
15 errors are known or determined. The following are written, where 
LOS; , LOS< , and LOS? are the x, y, and z vector components in the east, north, and 
up coordinate system of a line of sight to the ^th satellite: 

LOS* LOS'y LOSL 



H = 



LOS* LOS 



LOSz 



los; los" los; 
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y = 



A^-A(N!-N2)-c(dt,-dt 2 ) 
A^ 2 -A(N?-Nl)-c(dt,-dt 2 ) 

A^"-A(Nr-N5)-c(dt,-dt2). 



d = 



dx 

dy 

dz 



Vector d is a vector representing a baseline in an x, y, and z coordinate 
system, representing east, north, and up, respectively. Using the above equations, the 
following is also written: 

— ► -* 
y = Hd 

A solution to the above equation is written: 
d = (H r H) _, H r y , subject to the additional condition that: 



|| d ||= a distance between antennas 10 and 12. 

In another aspect of the present invention, double differencing, as 
known to those skilled in the art, is used to eliminate clock errors of receivers 16 and 
18. 

Locomotive 14 rides on a track having limited pitch. Maximum 
pitches range up to about 2 to 3 percent. Therefore, it is much easier to apply the 

constraint on || d || and to solve for integer ambiguities of N g p than in more general 

cases. Once d is determined, an attitude of locomotive 14 is determined at every 
epoch. By computing differences in the attitude, the attitude rate is computed. Clock 
errors and integer ambiguities are readily computed, so that a continuous 
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measurement of heading and heading rate is output from GPS receivers 16 and 18. In 

d , tan"' d 2 _ 

particular, heading is given by tan"'-^- , and pitch is given by ^ +rf2 • mus, a 

continuously updated attitude and attitude rate are available. It is only necessary that 
receivers 16 and 18 maintain phase lock with received signals from the GPS satellites 
20, 22... up to an nth satellite. If phase lock is not maintained cycle slip will occur at 
which time integer ambiguities will have to be recalculated. Microprocessor 24 
coupled to receivers 16 and 18 is used, in one embodiment, to compute at least one of 
an accurate heading and heading rate. As used herein, microprocessor 24 refers to 
controllers and processors, including microcontrollers, programmable logic 
controllers, input/output controllers, reduced instruction set circuits, application 
specific integrated circuits, logic circuits, and any other circuit, processor or 
microcomputer capable of processing the embodiments described herein. For 
example, results are displayed on display 26 or recorded on an internal or external 
device 28 such as a hard disk or magnetic tape or sent to a navigational computer . In 
one embodiment, device 28 is a transmitter configured to relay computed location 
and/or motion information to another location. 

A present absolute position of locomotive 14 measured in latitude and 
longitude is provided by GPS receivers 16 and 18. To resolve an initial integer 
ambiguity, determine clock offsets, or both, a database 30 is consulted to provide 
initial heading and grade of a track as a function of the latitude and longitude 
information. Given LOS values, three unknown parameters d x , d y , and d z are directly 

obtained by computation subject to the constraint that a magnitude of d is equal to a 
distance between antennas 10 and 12. 

In yet another embodiment and referring to Figure 2, track database 30 
contains measured or otherwise known track curvatures, indexed by location. A 
single GPS receiver and antenna, for example GPS receiver 16 and antenna 12, 
provides a position signal 54 (in addition to any other signals provided for other 
purposes). Using a location inferred from position signal 54, track database 30 is 
consulted to determine track curvature C at that location. 



WO 01/18558 



PCT/US00/24807 



Another aspect of the invention is detection of curves and reduction of 
track wear. Referring to Figures 3 and 4, when wheels 32 of locomotive 14 move 
over curve 34 on a track 36, resulting friction causes track 36 to wear down. A track 
is worn down, for example, when a crown (not shown) has worn off. By detecting 
curve 34 and dispensing a lubricant 38 to rails of track 36 while locomotive 14 
negotiates curve 34, wear of track 36 rails is reduced. The lubricant 38 may be 
dispensed from a lubricating system or any system for lubricating the rails of track 36, 
including lubricant dispenser 44. 

Referring to Figure 5, curves 34 are detected in one embodiment by 
utilizing a gyroscope 40 to detect an angular rotation rate co and a tachometer 42 to 
detect a train velocity v. Curvature C is then written C=colv. One or more lubricant 
dispensers 44 are configured to dispense lubricant 38 to track 36 rails. Specifically, 
lubricant dispensers 44 are operated by a microprocessor 24, for example, the 
microprocessor 24 of Figure 2. Microprocessor 24 is responsive to an angular change 
signal 46 related to co and a tachometer signal 48 related to velocity v to operate 
lubricant dispenser 44 to dispense lubricant 38 to track 36 rails when a predetermined 
magnitude of curvature C is exceeded. 

In another embodiment and referring to Figure 6, an angular change 
signal (not shown) is determined by microprocessor 24 utilizing signals supplied from 
a pair of GPS receivers 16 and 18 (such as those also shown in Figure 2) from signals 
received at antennas 10 and 12, respectively. This angular change signal is easily 
obtained because one attribute of attitude is a heading /i, and Ah/At=a>. Tachometer 
42 supplies a tachometer signal 48 in this embodiment, and microprocessor 24 
computes C=a/v. 

In yet another embodiment and referring to Figure 7, a lateral 
accelerometer 50 produces or generates a lateral acceleration signal a. Microprocessor 
24 is responsive to accelerometer 50, signal 52 and tachometer signal 48 to control 
lubricant dispensers 44. Lubricant 38 is dispensed when a predetermined magnitude 
of curvature C, as determined by a ratio a/v 2 , is exceeded. 
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Lubricant 38 is applied to reduce friction between wheels 32 and track 
36 rails when a magnitude of curvature C is exceeded by an amount to make track 
wear a concern. Accordingly, in one embodiment, lubricant dispensers 44 are 
configured to dispense lubricant 38 to wheels 32 rather than directly on track 36 rails. 

In embodiments in which measurement accuracy is unacceptably 
compromised by selective availability (SA), differential GPS (DGPS) is used. (SA 
results from a time-varying drift that either is, or that can be added deliberately to 
GPS satellite signals to limit absolute measurement accuracy for commercial and non- 
military applications.) 

Using one of receivers 16 or 18 in conjunction with microprocessor 24, 
another aspect of the present invention provides accurate measurements of distance 
traveled by the locomotive 14. However, it will be recognized that this aspect of the 
invention is also applicable to other land-based vehicles, as well. In one embodiment, 
a distance traveled is accurately determined even in the presence of SA with or 
without external DGPS connections. Samples of latitude and longitude from one of 
the receivers, for example, receiver 16, is sampled frequently enough so that errors in 
these samples are highly correlated. For example, samples are taken about every 
second. A distance between the samples is computed as a function of differential 
latitude and longitude. The difference Ad is written as follows: 

Ad = R[Alat 2 + cos 2 Qat) Along 2 f 2 

where Alat is a difference between latitudes of consecutive 

measurements; 

Along is a difference between longitudes of consecutive measurements; 

and 

R is the radius of the earth (about 3,440 nmi). 

A total distance traveled is a summation of Ad over successive 
measurements. This distance measurement is accurate even though absolute latitude 
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and longitude measurements are subject to large errors, because these errors are 
highly correlated between consecutive measurements and cancel out when differences 
are calculated. Even errors resulting from selective availability are readily cancelled 
by this method. 

From the preceding description of various embodiments of the present 
invention, it is evident that navigational parameters of a locomotive are efficiently and 
inexpensively measured utilizing two GPS receivers, and track wear is reduced 
through efficient use of lubricant. 

Even though a GPS satellite system is described, any satellite system 
or network may be utilized, for example, GPS, DGPS, and Glonass. 

Although the invention has been described and illustrated in detail, it is 
to be clearly understood that the same is intended by way of illustration and example 
only and is not to be taken by way of limitation. Accordingly the spirit and scope of 
the invention are to be limited only by the terms of the appended claims and their 
equivalents. 
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WHAT IS CLAIMED IS: 



1 . A method for determining at least one of motion and location 
parameters of a locomotive, said method comprising the steps of: 

determining a set of phase differences between satellite reference 
signals received by satellite receivers; and 

determining at least one of an accurate heading, heading rate, attitude, 
and attitude rate of the locomotive using the set of phase differences between the 
satellite reference signals. 

2. A method according to Claim 1 further comprising the step of 
determining d . 

3. A method according to Claim 2 further comprising the step of 
determining d as d = (H r H)"'H y , where: 

LOSl LOSl LOSz 



H = 



y = 



LOS* LOS y LOSz 

LOS" LOS y LOSS 

A^'-A(N!-N2)-c(dti-dt2) 
A^ 2 -A(N?-N')-c(dt,-dt2) 

A^ n -A(Nr-N5)-c(dt,-dt2). 



; and 



d = 



d* 

dy 

L dz. 
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4. A method according to Claim 3 wherein said step of 
determining at least one of an accurate heading, heading rate, attitude, and attitude 
rate of the locomotive further comprises the step of determining an attitude and an 

attitude rate of a locomotive using d , the heading using tan" 1 ^f- , and heading rate 



tan" 1 d 2 

us,ng ^F' 

5. A method according to Claim 1 further comprising the step of 
determining a track curvature, C. 

6. A method according to Claim 5 wherein determining a track 
curvature comprises the step of detecting an angular rotation rate co and a velocity v of 
the locomotive, wherein C=o)/ v. 

7. A method according to Claim 6 wherein detecting an angular 
rotation rate co and a velocity v of the locomotive comprises the step of detecting an 
angular rotation rate co using a gyroscope and a velocity v of the locomotive using a 
tachometer. 

8. A method according to Claim 6 wherein detecting an angular 
rotation rate co and a velocity v of the locomotive comprises the step of detecting an 
angular rotation rate co using received satellite signals and velocity v of the 
locomotive using a tachometer. 

9. A method according to Claim 5 wherein determining a track 
curvature comprises the step of determining a lateral acceleration a and a velocity v of 
the locomotive, wherein C= a I v 2 . 

10. A method according to Claim 5 further comprising the step of 
controlling dispensing of a track lubricant in accordance with the determined 
curvature value C. 
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11. A method according to Claim 10 wherein said step of 
controlling dispensing of a track lubricant further comprises the step of dispensing the 
lubricant when C exceeds a predetermined magnitude. 

12. A method according to Claim 1 further comprising the steps of: 
determining a position of the locomotive; and 

accessing a database of track heading and grade to determine a present 
track heading and grade at the determined position of the locomotive. 

13. A method according to Claim 1 further comprising the steps of: 
sampling latitude and longitude from the satellite receivers; and 
determining a distance traveled by the locomotive. 

14. A method according to Claim 13 wherein said step of sampling 
latitude and longitude from the satellite receivers further comprises the steps of: 

sampling where the distance between the samples is determined as 
Ad = R[Alat 2 + cos 2 fat) Along 2 ] 1 ' 2 

where Alat is a difference between latitudes of consecutive 

measurements; 

Along is a difference between longitudes of consecutive 

measurements; and 

R is the radius of the earth (about 3,440 nmi); and 

said step of determining a distance traveled by the locomotive further 
comprises the step of summing Ad over successive measurements. 

15. An apparatus for determining at least one of motion and 
location parameters of a locomotive to detect curves and reduce track wear, said 
apparatus comprising: 

12 
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at least two phase-locking satellite receivers configured to reference 
signals received from a set of satellites; and 

a processor configured to determine a set of phase differences between 
the reference signals received by said satellite receivers and at least one of an accurate 
heading, heading rate, attitude, and attitude rate of the locomotive using the set of 
phase differences between the reference signals. 

16. An apparatus according to Claim 15 wherein said processor 

-+ 

further configured to determine d . 

17. An apparatus according to Claim 16 wherein said processor 

_» -» -» 
further configured to determine d as d = (H r H)~'H y , where: 

"LOSL LOS' y LOSi" 
LOS, LOS 2 y LOSr 

9 

• • • 

LOS? LOS? LOS" 

A^'-A(N|-N2)-c(dti-dt2) 
A/-A(N?-N 2 )-c(dt,-dt 2 ) 

A^ n -A(Nr-N;)-c(dt,-dt 2 ). 

d," 
d y • 
_dz_ 
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18. An apparatus according to Claim 17 wherein said processor 
further configured to determine an attitude and an attitude rate of the locomotive 

i d 3 , fl ^ tan~ ! <f 2 

using d , the heading using tan" 1 -f- , and the heading rate using -p- = . 

d y ^d x +d y 

19. An apparatus according to Claim 15 wherein said processor 
further configured to determine a track curvature, C. 

20. An apparatus according to Claim 19 wherein to determine a 
track curvature, said processor configured to detect an angular rotation rate co and 
velocity v of the locomotive. 

21. An apparatus according to Claim 20 wherein to detect an 
angular rotation rate co and velocity v of the locomotive, said processor configured 
detect an angular rotation rate co using received satellite signals and velocity v of the 
locomotive using a tachometer. 

22. An apparatus according to Claim 20 wherein to detect an 
angular rotation rate co and velocity v of the locomotive, said processor configured 
detect an angular rotation rate co using a gyroscope and velocity v of the locomotive 
using a tachometer. 

23. An apparatus according to Claim 19 wherein to determine a 
track curvature, said processor further configured to: 

determine a lateral acceleration a and a velocity v of the locomotive; 

and 

determine track curvature C as: C = a / v 2 . 

24. An apparatus according to Claim 19 further comprising a 
device for dispensing a lubricant to a track. 
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25. An apparatus according to Claim 24 wherein said processor 
further configured to control said device dispensing the lubricant in accordance with 
the determined curvature value C. 

26. An apparatus according to Claim 25 wherein said processor 
further configured to dispense the lubricant when C exceeds a predetermined 
magnitude. 

27. An apparatus according to Claim 24 wherein said processor 
further configured to: 

determine a position of the locomotive; 

access a database of track heading and grade to determine a present 
track heading and grade at the determined position of the locomotive; and 

control said device dispensing lubricant in accordance with a curvature 
value C contained within the track database. 

28. An apparatus according to Claim 15 wherein said processor 
further configured to: 

sample latitude and longitude from the GPS receivers; and 
determine a distance traveled by the locomotive. 

29. An apparatus according to Claim 28 wherein said processor 
configured to determine a distance between samples as: 

Ad = R[Alat 2 + cos (l<*0 Along 2 ] 1 ' 2 . 
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